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ABSTRACT 

We have investigated the relationship between the kinematics and mass of young (< 
3 x 10 8 years) white dwarfs usi ng proper motions. O ur sample is taken from the 
colour selected catalogues of SPSS (lEisenstein et alj|2006l) and the Palomar-Green Survey 
dLiebert. Bergeron & Hol berg 2005), both of which have spectroscopic temperature and grav- 
ity determinations. We find that the dispersion decreases with increasing white dwarf mass. 
This can be explained as a result of less scattering by objects in the Galactic disk during the 
shorter lifetime of their more massive progenitors. A direct result of this is that white dwarfs 
with high mass have a reduced scale height, and hence their local density is enhanced over 
their less massive counterparts. In addition, we have investigated whether the kinematics of 
the highest mass white dwarfs (> 0.95 M© ) are consistent with the expected relative contri- 
butions of single star evolution and mergers. We find that the kinematics are consistent with 
the majority of high-mass white dwarfs being formed through single star evolution. 
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1 INTRODUCTION 

Despite the significant work on both the kinematics and mass dis- 
tribution of white dwarfs, very little work has addressed their con- 
nection. 

The kinematics of galactic white dwarfs have been studied on 
numerous occasions with several motivations. They have proven 
useful in attempts to unravel the evolutionary history and pro- 
genitors of the various classes of white dwarfs dSion et afl|l988t 
lAnselowitz et aljl999t) . Interest in white dwarf kinematics was also 
prompted by the suggestion that halo white dw arfs could provide a 
significant contrib ution to Galactic dark matter dOppenheimer et aT] 
1200 it iReidl l2005h . This effort has concentrated on the identifi- 
cation of halo white dwarfs and estimating the resultant density, 
which now appears to be a small contribution to the Galactic 
dark matter budget dPauli et aT1 l2006), Moreover, the mass distri- 
bution of the most common hydrogen rich (DA) white dwarfs has 
also been extensively investigated, particularly for white dwarfs 
with T > 10,000 K which are hot enough for their masses 
to be deduced spectroscopically from fits to their B aimer lines 
dLiebert. Bergeron & Holberd 120051: iKepler et all 120071; IVennesI 
1 19991) . The mass distribution shows a peak at O.6M0 due to the 
relative abundance of their lower mass progenitors with a tail ex- 
tending to higher masses formed from more massive progenitors. 

The connection between the galactic kinematics of a group of 
thin disk objects and their p rogenitors is largely due to the proces s 
of kinematic disk 'heating' dWielen|[l977l: iNordstrom et al.ll2004l) . 
The hot white dwarfs with short cooling ages we observe in the 
galactic neighbourhood today are formed from a wide range of pro- 
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genitor masses (~ 0.8-8 Mq) and hence have a wide range in age. 
We therefore expect high-mass disk white dwarfs to have a low 
velocity dispersion in comparison to low-mass disk white dwarfs 
w hose progenitors formed earlier. Thi s conne ction was suggested 
in iGuseinov. Novruzova & Rustamovl d 19831) who performed an 
analysis suggesting that white dwarfs with larger m asses have 
smalle r dispersions, however this was reinvestigated by ISion et al] 
dl988l) with a larger sample of 78 DA white dwarfs where no ev- 
idence for any correlation was found. This paper readdresses the 
connection between mass and kinematics with a greatly increased 
sample size. 

The outline of the paper is as follows: In sectionf2]we discuss 
the sample selection and the calculation of distances and proper 
motions. In section[3]we discuss how we estimate the kinematics of 
the sample without radial velocity in formation. We use two meth- 
ods, that of iDehnen & Binnevl J 1998h (section [3~TT l, and a Markov 
Chain Monte Carlo (MCMC) where we marginalise over the un- 
known radial velocity (section[372}. In section|4]we analyse whether 
the kinematics are consistent with single star evolution (SSE) both 
via analytic methods (section l4~H and simulations (section l4~2t . In 
section [5] we analyse whether the highest mass white dwarfs are 
largely formed through single star evolution or are the product of 
the merger of two lower mass white dwarfs. Finally, we discuss the 
implications of our findings on the scale height of white dwarfs in 
section[6] 

For the reader in a hurry, the primary result of this paper, the 
relationship between the mass of young white dwarfs and their 
velocity dispersion, is shown in figure [3] and discussed in section 
[3] The implied scale heights, the second key result, are then dis- 
cussed in section[6] These results have been checked using a Monte 
Carlo simulation of the formation and observation of an ensemble 
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of white dwarfs, which is described by flowcharts in figures [6}|8] 
in figure [6] the process of choosing stars is described, in figure [7] 
the process of placing them in the disk is described, and in figure[8] 
the process of determining the observability of the simulated white 
dwarf is described. 



2 SAMPLE 

We investigate only hydrogen atmosphere (DA) white dwarfs due 
to the relative simplicity of their spectra and the resultant security 
of the spectroscopic masses. The sample of DA white dwarfs is 
take n from two sources, the Palomar-Gr een (PG) white dwarf sur- 
vey dLiebert. Bergeron & Holberdl2005l) and the SDSS DR4 white 
dwarf survey ( Eisenstein et al. 20061) . The SDSS sample is much 
larger the the PG sample. The PG sample is included as a demon- 
stration that the results are secure, and not a result of systematics in 
SDSS, such as the complex selection of targets. For clarity we first 
discuss which types of white dwarfs we select, then discuss how 
the SDSS survey is dealt with, and finally how the PG survey was 
dealt with. The sample and its selection is summarised in table Q] 

Selected White Dwarfs: Both PG and SDSS are colour 
selected, eliminating the kinematic biases inherent in proper motion 
based surveys, and contain spectroscopic determinations of surface 
gravity, logg, and effective temperature, T e ff, obtained by fitting 
the profile of the Balmer lines. We restrict the sample to objects 
whose fitted T e g was between 13,000 K and 40,000 K, since log g 
appears to be systematically overestimated at low temperatures an d 
T e ff overestimated at higher temperatures (lEisenstein et al.lEood) . 

The fitted log g and T e g are converted to masses and ages us- 
ing the models of the carbon core wh ite dwarf cooling mo dels of 
iFontaine. Brassard & Bergeronl (l200lh below 30,000 K and lWoodl 
19951) with thick hydrogen layers of fractional mass 10 4 above 
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30,000 K \ White dwarfs with inferred masses less than 0.47 Mq 
are instead assumed to have heliu m cores whose masse s and ages 
are calculated from the models of ISerenelli et al.l d200lh . Only ob- 
jects with cooling ages below 3 x 10 8 years are included in the 
sample to avoid significant kinematic heating after white dwarf for- 
mation. The requirements of cooling age below 3 x 10 8 years and 
T cff above 13,000 K are competing. Above 0.60 M the WDs cool 
more slowly and thus the age limit is used, while below 0.60 Mq 
the temperature limit is used. 

White dwarfs previously discussed in the literature as known 
members of binaries w ere removed from the sampl es. 

SDSS SURVEY dElSENSTEIN ET AlJ|2006|) : Many of the 
SDSS spectra have low signal-to-noise ratios and hence large errors 
on their fitted log g and T c ff . To ensure accurate masses and pho- 
tometric distances only objects whose spectra had a signal-to-noise 
ratio larger than 10 are included. The grid of model atmospheres 
fitted in the SDSS catalog extends only to log g = 9, and thus, for 
object s at this limit, the refitted log g and T c s given in lKepler et alj 
d2007h were used. 

Photometric distances to the white dwarfs in SDSS are calcu- 
lated by minimising 

X 2 = (m i -[M i (logg,T cH ) + 

i—(u,g,r,i,z) 

A g a, + 5 log d-5]) 2 /a 2 (1) 

where mi and at are the 5 band SDSS photometry and 
their errors, Mi are the model absolute magnitudes, A g a,i is 
the reddening and d the distance in parsecs. The photomet- 
ric <Tj is the quoted photometric error in SDSS each band 
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Figure 1. \ 2 per degree of freedom (DOF) for the fitted photometric dis- 
tance of the 1443 SDSS DA white dwarfs considered. A x 2 function with 
3 DOF is plotted as the dotted line. Beyond \ 2 = 5 the white dwarfs are 
rejected. 



added in quadrature to a system atic error of (u,g,r, i,z) = 
(0.015,0.007,0.007,0.007,0.01) dKleinman et "all 1 20041) . Model 
absolute magnitudes are taken from the atmospheric models pro- 
vided by BergerorQ. A g a,i is the product of Ry =3.1 extinction in 
each band of (a u , a g , a r , a it a z ) — (1.36,1.00,0.73,0.55,0.39) 
and the overall extinction A g , which is constrained to lie 
between zero and the value of galactic extinction map of 
ISchlegel. Finkbeiner & Davisl dl998l) at the position of the object 
considered. 

The resulting distribution of \ values calculated by minimis- 
ing equationQ]is plotted in figure[T] It closely resembles a \ 2 dis- 
tribution, but with an extended tail. Objects with reduced \ 2 larger 
than 5 were removed from the sample, most of these objects show 
an excess towards the redder photometric bands, indicating they are 
in binaries with a cooler white dwarf companion. Errors in the pho- 
tometric distance are taken to be the A^ 2 = 1 surface added in 
quadrature to the distance errors introduced though the uncertainty 
in log g and T e ff . 

Pro per motions f or the SDSS sample are taken from the cata- 
logue of iMunn et al. (2008). These proper motions are calculated 
from the USNO-B1.0 plate positions re-calibrated using nearby 
galaxies together with the SDSS position so that the proper mo- 
tions are more a ccurate and absolute . By measuring the proper mo- 
tions of quasars IMunn et al] d2004l) estimates that the la error is 
5.6 mas yr -1 . 

PG SURVEY: For 132 stars in the PG survey, SDSS photome- 
try was available and the same method was used as for SDSS stars. 
For the remaining objects the PG c atalog photometric distances 
were u sed. These were estimated in iLiebert. Bergeron & Holberd 
d2005l) from comparison of the V b and magnitude with t he pre - 
dicted My from the same models of iHolberg & Bergeronl d2006l) . 
Comparison of the stellar distances given by the two methods gives 
a standard deviation of 7 per cent. The majority of this error is ex- 



1 Available from |http : / / www .astro . umontreal ■ ca/ -bergeron/Cooling Model 
uses r e sults from from|H olberg & Bergeron (2006). Kowalski & Saumon 
1 2006 ). [Tremblav, BergeroiT^Giaiuiinasni'oill) and lBergerorTetal 



© 201 1 RAS, MNRAS 000,[TJ[T4] 



White Dwarf Kinematics vs Mass 3 



PG SDSS 



Number of DA White Dwarf s with 


299 


6926 


good photometry not known to be binaries 


of these number with signal-to-noise > 10 


299 


3125 


of these number with 13, 000A" < T off < 40, 000K 


215 


1555 


of these number with age < 3 X 10 s yrs 


21 1 


1491 



Distance source: 






Liebert. Bergeron & Holbere (2005) 


79 





SDSS Photometry 


132 


1491 


of these number rejected with \ 2 > 5 





48 


Proper Motion Source: 






Munn et al. (2008) 


153 


1443 


PPMXL 


54 





Manual measurement from POSS I/II 


4 






Table 1. Summary of sample 



pected to be in the PG survey distances and hence a conservative 
10 per cent error was applied to these. 

Proper motions for PG white dw arfs that ap p ear i n 
SDSS are taken from the catalog of iMunn et al.l 
For the remaining objects, the PPMXL proper motion 
was used w here available, which has typical la error of 
~8 masyr -1 teoeser. Demleitner & S chilbach 2010h . 

Finally 4 objects in the PG sample have no reliable PPMXL 
proper motion, primarily due to a spurious matching of objects be- 
tween epochs. For these, the proper motion was calculated directly 
between the scanned POSS-I and POSS-II plates. The proper mo- 
tion was measured relative to nearby faint stars of similar magni- 
tude corrected for galactic rotation (see Section l3TTl . Typical errors 
estimated from the proper motions of stars of similar magnitude to 
be 11 masyr -1 . We emphasize that only 4 of 1491 white dwarfs 
use this method, and none have mass above 0.95 Mq analyzed in 
more detail in Section|5] 

FINAL SAMPLE: The resulting sample of 1443 SDSS and 211 
PG white dwarfs contains young DA white dwarfs with reliable 
masses, proper motions and photometric distances. The mass distri- 
bution of the samples is shown in figure|2] The process of construct- 
ing the sample together with numbers of objects is summarised in 
tableU 



3 KINEMATICS WITHOUT RADIAL VELOCITIES 

We now turn to calculating the mean velocity and the velocity dis- 
persion for our sample. While radial velocities are required to com- 
pletely determine the kinematics of an individual object, bulk kine- 
matic properties such as the mean velocity and the velocity disper- 
sion can be determined from only transverse motions. 

We use two methods to do so, the frequentist method used in 
section [3~71 and a Markov Chain Monte Carlo in section [3~2l Both 
methods give similar results which are summarised in table [2] 



3.1 Method of Dehnen and Binney (1998) 

The method used here is adapted from iDehnen & Binnevl d 19981) . 



First the observed proper motions in galactic coordinates. 
Hi, bB , are corrected for Galactic rotation through 

Acos(2f) - B 



obs 



fli 
fit, 



obs 
IJ'i 

/if 8 + Asm{2£) cos b sin b. 



and 



(2) 
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Figure 2. Mass distribution of the samples of SDSS (black line) and PG 
(dashed red line) white dwarfs after the cuts described in the text. Inset 
graph shows the same data on logarithmic axes. 

using A = 14.82 kms" 1 kpc" 1 and B = -12.37 kms -1 kpc" 1 
dFeast & Whitelocklll997h . In galactic coordinates where the com- 
ponents are directed towards the galactic centre, in the direction of 
galactic rotation, and towards the north Galactic pole we observe 
the velocity 



4.74d 



kms 



—He sini cos b — fib cos £ sin b 
Hi cos i cos b — [ib sin i sin b 
lib cos b 

(3) 

with d in kpc and proper motions in masyr -1 . This is the pro- 
jection of the velocity V onto the sky plane though the projection 
matrix 



V± = A ■ V , A = I-i 

where r is the unit vector to the star. 

Next the quantity S 2 is formed through 



S 2 (Vo) = (\V ± - A- Vb| 



(4) 



(5) 



Under the assumption that the positions of the observed objects are 
uncorrected with the velocity, then the choice of Vb that minimises 
S 2 is the mean velocity. Also S 2 at the minimum is a measure of 
the dispersion of the group. 

IDehnen & Binnevl dl998h then calculate all independent six el- 
ements of the dispersion tensor. Unfortunately, this entails estimat- 
ing nine parameters which limits its use to samples with large num- 
bers of objects. This would require excessively wide bins for the 
high-mass region where there are few objects. Instead we choose 
to make further assumptions about the objects' velocities in or- 
der to reduce the number of fitted parameters. The mean velocity 
of each group towards the galactic centre and the north Galactic 
pole is simply a result of the solar motion and we take these to be 
10.00 kms -1 and 7.17 kms -1 respectively jDehnen & Binnevl 
1 19981) . The mean velocity in the direction of galactic rotation, Vb, 
is kept as a free parameter since in addition to the solar motion this 
varies between groups due to asymmetric drift. We also assume that 
the dispersion tensor takes the form 



a — aidiag I 1 



J_ J_ 

1.4' 2.2 



(6) 



which is accurate for main sequence stars in the solar neighbour- 
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hood dDehnen & Binnevl 19981) . This reduces the number of param- 
eters for each group to the asymmetric drift Vo and the normalisa- 
tion of the dispersion tensor ai . 

Vo is calculated by minimising equation [5] and then ai is es- 
timated though a Monte-Carlo simulation: Since S 2 is a measure 
of the dispersion, an initial estimate of a'{ is taken to be S 2 , and 
a set of simulations is performed where a new velocity is chosen 
for each white dwarf at its position in the sky from the isothermal 
distribution with the assumed dispersion tensor and the calculated 
mean velocity. The error in tangential velocity, assumed to be Gaus- 
sian, is added to this. The set of simulations produces a distribution 
of S 2 values, and a\ is iterated until the mean S 2 corresponds to 
the value calculated from observations. S 2 is almost proportional 
to a\ when errors in tangential velocity are neglected and so the 
error in a 2 is estimated from the distribution of S 2 scaled by this 
proportionality constant. 



3.2 MCMC Estimate 

In addition, a Markov Chain Monte Carlo (MCMC) likelihood 
based estimate of the kinematic parameters was obtained. We use 
uninformative fiat priors for the fitted parameters. 

We denote the probability that the velocity of the zth object 
was V to be P(V\Di, o - ,) where Di — (l,b, d, /J,b) is the 
data for the ith object together with the corresponding errors ai. 
Hi and are the values corrected for galactic rotation by equation 
[2] Under the assumption that positions are uncorrelated with ve- 
locity then the distribution function is a function only of velocity: 
f(V). In addition, in what follows we do not consider the posi- 
tions, but instead focus on the kinematics through the velocity V. 
Under these assumptions the overall likelihood for a set of obser- 
vations of a group of white dwarfs is 

c = n[dVf(y)P{v\D i ,* i ) (7) 

=► log£ = ^log f dVf{V)P(V\Di,<n) (8) 

i 
i 

In calculating the likelihoods, d, we assume a Schwarzschild dis- 
tribution function, and normally distributed error in proper motion. 
The unknown radial velocity is integrated over analytically. Ex- 
plicit expressions for d are given in appendix lAl 

Again, the dispersion tensor and mean were constrained to re- 
duce the number of parameters. We use flat priors on the dispersion 
and asymmetric drift. The expression for the likelihood was used 
to calculate the maximum likelihood estimate of the dispersion ten- 
sor, while errors were estimated from a MCMC using Metropolis- 
Hastings sampling. When the constraints on the dispersion tensor 
and mean velocity were relaxed this did not substantially alter the 
results, aside from the larger errors, particularly in the underpopu- 
lated bins due to the reduced degrees of freedom. In particular, the 
results are insensitive to allowing vertex deviation. 

The fitting results for the SDSS and PG samples are sum- 
marised in table [2] and plotted in figure [3] In addition, in figure 
[4] the raw transverse velocities measured from the proper motions 
for three groups of white dwarfs are shown. The lowest mass white 
dwarfs, AI < O.45M0, are expected to be predominantly formed 
through binary evolution and have a binary white dwarf partner. 
This potentially introduces errors into their photometric distances 
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Figure 3. Dispersion of SDSS (black) and PG (red) white dwarfs against 
mass calculated using the using the method described in section [3~2l Each 
bin is plotted at its mean mass. 
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Figure 4. Cumulative transverse velocity distribution of the combined 
SDSS and PG surveys. Low mass white dwarfs (0.5 Mq ^ M < 
0.75 Mq, with M = Mi + M2) as solid black, high mass white dwarfs 
(M > 0.95 Mq) as dashed red, and intermediate mass white dwarfs 
(0.75 Mq < M < 0.95 M©) as dotted green. 



and so we do not consider them beyond simply stating the fitting 
results in table [2] 
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Table 2. Kinematic fitting results from the PG and SDSS samples described in section [2] using the methods of sections ^. 1 | and [T2l Mi ow and A/high are m 
units of Mq, while a\ and V are in km s~ . N is the number of white dwarfs in each mass bin. 



PG SDSS 



Mow 


A/high 


JV 


Dehnen & 


: Binnev (19981 


MCMC 


A ; 


Dehnen & Binnev (1998) 


MCMC 










V 


<y\ 


V' 




ci 


V 


Cl 


V 


0.30 


0.40 


5 


47 ± 12 


22 ± 13 


48 ± 10 


18 ± 12 


70 


53 ±3 


33 ±4 


40 ±3 


34 ±4 


0.40 


0.47 


20 


49 ±6 


27 ± 7 


49 ±7 


28 ±6 


62 


68 ±4 


38 ±6 


70 ±6 


38 ± 6 


0.47 


0.55 


35 


47 ±4 


18 ±6 


51 ± 1 


17 ±4 


333 


56 ±1 


34 ±2 


57 ±2 


34 ±2 


0.55 


0.60 


53 


37 ±2 


17 ±3 


40 ±3 


18 ±3 


482 


46 ±1 


20 ±1 


45 ±1 


21 ±1 


0.60 


0.65 


51 


37 ± 2 


16 ±3 


34 ±2 


15 ±2 


239 


33 ±1 


20 ±1 


31 ±1 


20 ±1 


0.65 


0.75 


23 


33 ±3 


15 ±4 


34 ±4 


14 ±5 


91 


26 ±1 


16 ±2 


28 ±1 


15 ±1 


0.75 


0.85 


9 


16 ±2 


11±4 


17 ±3 


11 ±4 


30 


16 ±1 


11 ±2 


19 ±2 


11 ±2 


0.85 


0.95 


10 


12 ±2 


15 ±2 


12 ±2 


13 ±2 


28 


18 ±1 


12 ±2 


19 ±2 


11 ±2 


0.95 


1.44 


5 


22 ±5 


14 ±7 


24 ±6 


12 ±6 


9 


19 ±3 


9±5 


24 ±5 


9 ± 6 



4 EXPECTATIONS FROM SINGLE STAR EVOLUTION 



4.1 Analytic 

In this section we describe the reasons for the relationship WD 
mass and dispersion within a simple analytic model, before moving 
onto the more complex Monte Carlo simulations of section l4~2l 

Within the framework of single star evolution (SSE) an en- 
semble of white dwarfs with the same mass would be expected 
to have a dispersion cr(tTOT), where a(t) is the disk heating re- 
lation, and tTOT is the total age of the white dwarf including its 
precursor lifetime (i.e. total pre-white dwarf stellar lifetime). Here 
iTOT will be given by £tot = (wd + isSE(Mi(M W D)) where 
two is the cooling age of the white dwarf and tssE(Mi(MwD)) is 
the total precursor lifetime, which is a function of the white dwarf 
mass through the initial-final mass relation (IFMR) Mj(M/). Two 
components of this prediction are particularly uncertain: the disk 
heating relation and the IFMR. We discuss these now. 

The best constraints on the IFMR come from open clusters. 
Spectroscopic fits of the masses of white dwarfs give the final mass. 
The initial mass is estimated using isochrone fitting to the main 
sequence turnoff to calculate the age of the cluster, which finally 
allows the corre sponding initial mas s to be inferred using the pre- 
cursor lifetime dCatalan et al .1120081) . This method has succeeded 
in producing IFMRs with a typical uncertainty of less than 20%. 
The strong dependance of the precursor lifetime on mass however 
makes this a considerable uncertainty in the dispersion relation. 

The most accurate data on the disk heating relation is given in 
Nordstrom et al. I d2004l) from an analysis of F and G dwarfs with 
radial velocities and Hippar cos data, although this data still permits 
a range of heating models ( Seab roke & Gilmorell2007l) . However, 
fo r consistency, we inste ad use the disk heating models estimated 
in ljust & Jahreifil feOlOh . since we also use their companion star 
formation histories. 

The effect of these model uncertainties are shown in figure [5] 
for the models described in table|3] Qualitatively the results appear 
to agree with the predicted relations: for white dwarfs more massive 
than 0.75 Mq the white dwarf progenitors precusor lifetime is short 
and there is little dependance of the kinematics on mass. Below 
0.75 Mq the dispersion sharply increases as the progenitor lifetime 
approached 1 Gyr and longer where the disk heating is significant. 

However, while qualitatively the results in figure[5]are consis- 
tent, there is quantitative disagreement. To assess this disagreement 
we turn to a more sophisticated Monte Carlo treatment. 




10 r 



OS; i i i i i i I 

0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

M/M 

Figure 5. Predicted dispersions from models A-D, as dashed blue, solid 
blue, dashed green and solid green respectively. Also plotted are SDSS 
white dwarfs (black error bars) with the data taken from the MCMC col- 
umn of tablefS] 



4.2 Monte Carlo 

As a quantitative check of our results in section |4~T1 we have per- 
formed a Monte Carlo simulation of the production, kinematics, 
and observation of the white dwarfs in the solar neighbourhood, 
as described in this section. We also describe the simulated selec- 
tion and observation of these white dwarfs by SDSS and PG. We 
perform this simulation to assuage fears that our results could be 
impacted by effects such as selection biases. 

This process is somewhat involved, and so for clarity it is sum- 
marised in the flow charts in figures [6}{8] The final results of the 
Monte Carlo simulation are compared with the white dwarf sample 
in figure [Tol 

PICKING Stars: The initial mass was drawn from a Kroupa 
IMF and one of two star formation histories (Table [3}. If this re- 
sulted in a white dwarf at the present time with an age less than 
3 x 10 8 years, and a temp erature betwe en 13,000K and 40,000K 
using the cooling models of lWoodl d 19951) as explained in sectionfJ] 
then it was included in the simulation. See figure|6]for synopsis. 

PLACING Stars IN Disk: If a star has been included in the 
simulation, it is given a velocity dispersion taken from the previ- 
ously described disk heating models of tableland axis ratios of the 
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Table 3. Model input parameters for the models of single star evolution (SSE). 



Model a(t) / kms" 1 Mi(M WT >) / M *ssb(A^) / Gyr SFR(i^_ 

A 66 ( 'ot+i2 yr ) 1/2 .l From Hurley. Pols & Tout (2000) . solar metallicitv. 3.25 El 

B 62 ( 'o 2 32Yio yr ) 1/2 - From Hurley. Pols & Tout (2000) . solar metallicitv. 7.68 cxp(-f./8 Gyr) G3 



0.32+10 ) 

C 66 ( °'ot+i2 yr ) 1/2 ^ From Catalan et al. (2008) From Girardi et al. (2000^ 3.25 El 

D 62 ( "'o^t+S^ ) 172 ® From Catalan et al. (2008) From Girardi et al. (2000'P 7.68 exp(-f./8 Gyr) 



In units of Mp) pc 2 Gyr . Not used in the an alytic SSE simulatio n of section l4~T 



Just & Jahreii 



Just & J ahreifi (2010) model C. Disk age 12 Gvr. lGirardi et al. 



201C ) model D. Disk age 10 Gvr.lGirardi et al 
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Figure 6. Flowchart illustrating the process of simulating white dwarfs 
formed from single star evolution (SSE). If a star reaches the final stage, 
then it is placed in the disk using a process described by the flowchart shown 
in figure[7] 



velocity ellipsoid of 1:1/1.4:1/2.2 jDehnen & Binnevl 19981) . Its ve- 
locity in the disk was drawn from a Gaussian with these widths and 
it was placed in the plane of the Galaxy using a radial exponential 
disk with a scale length of 2.5 kpc. Since the furthest > O.47M0 
WD projected into plane is less than 1 kpc, only WDs placed within 
this distance are simulated further. 

For an isothermal population the vertical position, z, and ve- 
locity, Vz , are given by 

oc exp(-i ) 2 /2^)exp(-$ 2 (2)/a 2 ) , (10) 

where $ z is the gravitational potential. Each star's velocity is thus 
drawn from a Gaussian with standard deviation given by the previ- 
ously calculated a z , while z is chosen by first drawing <J> 2 (z) from 



an exponential distribution with scale a t, and then inverting this to 
calculate z. We use the mass models of lHolmberg & FlvnrJ d200oT) 

for $z(»- 

This process of placing white dwarfs in the local galactic disk 
is summarised in figure [7] 

White Dwarf Observability: As a result of this process, 
each white dwarf has an assigned galactic position and velocity, 
together with its mass and age. It is then assessed whether it is 
likely to be observed in either the SDSS or PG survey as follows: 
First its galactic position is translated to a right ascension, RA, and 
declination, dec, and, unless this falls on one of the PG plates or 
the SDSS DR4 spectroscopic plates, the probability of observation 
is zero. 

For white dwarfs in the PG surve y the apparent U and B mag- 
nitude is calculated from the models of lHolberg & Bergeronl ( l20fJq) 
with a 0.27 mag error added to each to mimic the photometric er- 
rors in PG dLiebert. Bergeron & Holberg||2005l) . If it is bluer than 
U — B — —0.46 and brighter than the B band magnitude limit for 
the PG plate on which it lies then it is considered observed. 

For SDSS the spectroscopic targeting is more complex 
jKleinmanetall 120041) , and the strategy was to construct an em- 
pirical observational probability for a star at each magnitude and 
colour. A four dimensional table of probability of spectroscopic 
follow up was constructed in (r, u — g, g — r, r — i) grouped in 
0.2 mag bins from the SDSS DR4 clean photometry. The expected 
spectroscopic signal-to-noise was calculated using a quadratic least 
squares fit to the observed signal-to-noise ratio as a function of 
g-band magnitude together with normally distributed scatter in 
signal-to-noise with standard deviation of 1.7. If the signal-to-noise 
ratio was greater than 10 it was included in the mock sample. 

Finally, measurement errors in mass of 0.03 Mq and proper 
motion errors of 5.6 mas yr~ are introduced. 

The process of assessing if each white dwarf is observed by 
the PG or SDSS surveys is summarised in figure [8] In all simula- 
tions we simulate a total of ~ 2 x 10 11 objects. 

MONTE Carlo RESULTS: The results of this simulation are 
shown in figure [10] As a further check that the simulated white 
dwarfs have the correct kinematics we plot the distributions in the 
U, V and W directions (directed towards the galactic centre, in the 
direction of galactic rotation, and towards the north Galactic pole 
respectively) in figure [9] 

The results of the singe star evolution (SSE) simulation, de- 
scribed in this section, closely agree with the observations, mod- 
ulo the normalisation factor. We do not concern ourselves with this 
overall normalisation, however the normalisation factor is typically 
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(a) SDSS. 



(b) Palomar-Green (PG). 



Figure 8. Flowchart illustrating the process of simulating whether white dwarfs are observed. This process is undertaken if a star has reached the final stage 
of the flowchart shown in figure[7] 



< 2. The simulation also does not produce white dwarfs below 
m 0.47 Mq, which are generally expected to form through binary 
evolution. As may be expected from the analytic models plotted in 
figure [5] the models in table [3] all produce white dwarfs that rea- 
sonably closely explain the observed samples and their kinematics 
and so we only plot the results of only one representative model in 
figure [TO] 



5 EXPECTATIONS FROM BINARY STAR EVOLUTION 

It has been suggested that the majority of high-mass white dwarfs 
were formed from merger s of binary white dwarfs, both on the basis 
of their number density l lLiebert. Bergeron & Ho lberg 20051) an d 
a possible peak at 1 Mq i Vennesl 19991) . To test this hypothesis 
we use two binary evolution codes (discussed in section I5.lt to 
perform binary population synthesis (described in section [5~2l , and 
ultimately what fraction of the sample is likely to have had a binary 
WD progenitor (section l5"3l . 



5.1 Binary Evolution Codes 

To address the considerable uncertainties in binary evolution, two 
binary evo lution codes were us e d. Spe cifically, the BSE code de- 
scribed in iHurlev. Pols & Toutl ||2000|) . and the SeBa code de- 
scribed in iNelemans et al. 1 2001 ). Both codes use the same ap- 
proach to modelling binary evolution: semi-analytic fits to the 
structure and evolution of isolated stars are combined with prescrip- 
tions for interactions between the stars. 

There are four key initial conditions that govern the evolution 
of a binary: the initial primary mass Mu, the initial secondary mass 
M2i (or equivalently the mass ratio (ft — Mn/D/hi), the initial 
binary semi-major axis ai and the initial eccentricity e t . 

One slice through the four-dimensional space of initial condi- 
tions (Mn, (ft, Hi, ei) showing those conditions which result in the 
merger of a pair of white dwarfs is shown in figure [TT1 

The differences between the BSE code and SeBa code in fig- 
ure[TT]are striking, and are largely due to the different binary evolu- 
tion prescriptions, and in particular the treatment of the Roche lobe 
overflow (RLOF) and common envelope (CE) phases. 

For the specifics of the treatment of the RLOF phase and 
its treatment in the BSE and SeBa codes we refer the reader 
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Figure 10. Comparison between the observed white dwarfs as the solid lines, and the Monte Carlo simulations of single star evolution described in section POl 
as the dashed lines. The upper panels show the PG survey and the lower panels the SDSS survey. The rightmost panel shows the cumulative transverse velocity 
distribution. In this panel colors are as in figure|4] Specifically low-mass white dwarfs (0.5 Mq Mi + M2 < 0.75 Mq) are plotted in black, high-mass 
(Mi + M2 > 0.95 Mq) white dwarfs in red, and intermediate mass white dwarfs (0.75 Mq ^ Mi + M2 < 0.95 Mq) in green. The simulation plotted is 
model D from table [5] 



Model 


Evolution Code 


CE Prescription 


«cbA 


7 


i 


BSE 


act 


2 




ii 


BSE 


a a 


1 




iii 


SeBa 


yet 


2 


1.5 


iv 


SeBa 


act 


2 





Table 4. Summary of the four b inary evolution models co nsidered. 
The BSE code is that des cribed by iHurlev Pols & TouJ feOOd) , and the 
SeBa code is described in Nelemans et al The common envelope 

(CE) prescription describes how the two phases of common envelope evo- 
lution are treated. For example 7a describes treatment of the first phase 
through the 7 parameterisati on and the second thro ugh the a parameteri- 
sation. We refer the reader to iNeleman s et al.l |200lj) for the definition and 
descriptions of these parameterisations. 



to lHurlev. Pols & Toutl feoOOl) and INelemans et all d200lh respec- 
tively. There is also considerable uncertainty in the treatment of 
the important CE evolution phase. The most fundamental differ- 
ence between the codes is the treatment of the first phase of mass 
transfer. BSE uses the most commonly used prescription for com- 
mon envelope evolution known as the a parameterisation, for both 
phases of mass transfer. SeBa howeve r, by default use s the 7 pa- 
rameterization. We refer the reader to INelemans et al.l (1200 lh for 
the details of these parameterisations. 

To assess the result of the considerable uncertainties in binary 
evolution on the merger time distribution, and therefore the resul- 
tant velocity distributions, we have used four models across the two 
binary evolution codes. These models are summarised in table|4] 



5.2 Binary Population Synthesis 

We now describe our method of binary population synthesis. 

We use the sam e dist ribution s in the par a meters 



(Mn, qi, a\_, e{) as iHanl { 19981) and INelemans et alj 
with the exception o f the IMF for which w e use a Kroupal 



2001) 



2001) 



IMF as opposed to a lMiller & Scald J19791) IMF. For reference the 
probability distributions are: 



P(Mii) oc My™ 
P(qi) oc const. 
P(log di) oc const. 
P(ei) oc a 



0.8 < Mi; «: 10 , 

0< 9 <1, (ID 

< log ail R Q ^ 5 , 
«: ei < 1 . 



Our approach to simulating the results of binary star evolution 
is to first produce a 4-dimensional grid of binary simulations in 
the parameters (Mu, qi, eti, ei). Grid points were linearly spaced in 
Mu between 0.8 and 10 Mq, linearly spaced in qi between and 
1, logarithmically spaced in on between 1 and 10 4 Rq, and linearly 
spaced in ef between and 1. The grid size used was a 25 x 25 x 
50 x 10 grid in (Mu, qi, <n, ej), respectively. With this choice of 
grid combined with the distributions in equation QT] the population 
synthesi s is particularly simple: an initial primary mass is drawn 
from the lKroupj J200 lb IMF and a random binary from the closest 
corresponding (qi, ai, ei) slice is chosen. In all simulations a total 
of ~ 10 13 objects are places in the disk. 

The process of simulating stars formed from binary evolution 
is summarised in figure [121 

In what follows we concern ourselves with the merger of 
CO+CO white dwarfs, since these are the mergers proposed to 
result in > IMq white dwarfs. Thus, in figure [HI we plot the 
rate at which pairs of white dwarfs with sub-Chandrasekhar total 
mass merge as calculated from our binary population synthesis of 
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R, from number density dv oc exp(— i?/ Ro)RdR 



Get random orientation in galactic plane from 
sun from uniform distribution 




Get WD from SSE/BSE with mass M WD 
and age since star formation <tot 



Get cr z from age-dispersion relation using *tot 



Get random v z and <& z (z) assuming isothermal 
distribution f z oc exp(— v^/2af) exp(— Q z (z)/a1) 



Calculate z by inverting the <E' z (z) given by 
the mass model of Holmbcrg & Flynn (2000) 



distribution with a from age-dispersion relation 



Add solar motion and asymmetric drift of (v y ) 
= -o|/80 km s^ 1 (Dehnen and Binney, 1998) 



T 



Is the white dwarf observable at this position? 
See observability flow chart 



Figure 7. Flowchart illustrating the process of placing white dwarfs in the 
galactic disk and picking their velocity. This process is undertaken if a star 
reaches the final stage of the flowchart shown in figure [7] If a star reaches 
the final stage of this flowchart, the observability is finally determined using 
the algorithm described in the flowchart shown in figure [8] 



the four models in table [4] Note that the overall normalisation can 
be very different. In particular, model ii uses a relatively efficient 
CE prescription with qceA = 1 for both phases of mass trans- 
fer. This in turn results in a smaller range of initial separations 
that will ultimately result in a gravitational radiation driven WD 
merger. Despite the differences in overall rate between the mod- 
els they all display a similar distribution of merger times. This is 
because, apart from at early times, the merger time is dominated 
by the time to merge by gravitational radiation. This is a strong 
function of separation, a, specifically £gw oc a 4 . As a result, at 
late times, the merging WDs originally formed a narrow range in 
separation at WD+WD birth. Approximating this as a power law, 
^cco £ leads to a merger rate # = oc t~ {3 - c)/i , and 

da D at da at 

so for a wide range of e the me rger rate declines as ^ ~ t^ 1 
JlVlaoz. Sharon & Gal- Yarr]|2010h . 
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Figure 9. Histograms showing the agreement between the observed and 
simulated velocity distribution in U, V, W directions of SDSS WDs. The 
black line is the observed distribution, while the dashed red line is the dis- 
tribution of the SSE simulation for model C. Zero radial velocity is artifi- 
cially assumed, and number of simulated WDs is normalised to the number 
observed. U is directed towards the galactic centre, V in the direction of 
galactic rotation, and W towards the north Galactic pole. 
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Figure 11. Comparison of the WD+WD merger outcomes from the 
SeBa and BSE codes with their default prescriptions for binary evolution. 
All simulations use an initial mass ratio of qi = 0.5 and eccentricity of 
e\ = 0. Green corresponds to CO+CO, red He+CO and blue He+He. The 
lighter green are sub-Chandrasekhar (Mi + M2 < 1.4 Mq) mergers, and 
the dark green super-Chandrasekhar. 
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Figure 12. Flowchart illustrating the process of simulating white dwarfs 
formed from binary star evolution. 

10 E i 




t/Gyr 

Figure 13. Merger rates of CO+CO WDs with sub-Chandrasekhar total 
mass following a burst of star formation. The error bars are purely statis- 
tical due to the finite size of the simulated binary grid. Lines are models i-iv 
as solid line, dotted line, dashed line, and dash-dot line respectively. The 
models are described in table 0] SNuM = 1/(100 yr) /(10 10 Mq ). 



5.3 Proportion of high-mass White Dwarfs Formed in 
Mergers 

To assess the possible proportion of high-mass white dwarfs 
that formed through mergers, the CO+CO merger products with 

0. 95 M Q ^ Mi + M 2 < 1.4 M from the binary population syn- 
thesis, are subjected to the same process as the single population 
synthesis results i.e. they are placed locally in the disk according 
to the method summarised in figure |7]and their observability in the 
SDSS and PG samples assessed according to figure[8] 

We assume that no mass is ejected during the merger so that 
resultant white dwarf has mass Mwd = Mi + M2 ■ We also assume 
that the merger reheats the white dwarf sufficiently that the white 
dwarf has a cooling age of 

two = iform — t 

where iform is the time prior to the present at which the binary ini- 
tially formed, and t me r g e is the length of time it took for the merger 
to occur, including the precursor lifetime. The resulting cumulative 
transverse velocity of 0.95 Mq s£ Mi + M 2 < 1.4 M CO+CO 
merger products are shown in figure H4] 

In figure [14] and the following we have combined the PG and 
SDSS samples to improve the statistics. We combine the Monte 
Carlo results by the empirical proportions of WDs in this sample 

1. e. the observed PG to SDSS ratio of 5:9. Note however there 
is a possible discrepancy between the two samples in this high 
mass bin. In particular the SDSS sample has few low velocity 
(< 14 kms -1 ) white dwarfs (see the bottom right panel of fig- 
ure [TO), and this results in a 12% probability that they are drawn 
from the same distribution. 

The distribution of transverse velocities in figure[T4lshows that 
despite the uncertainties in binary evolution resulting in very differ- 
ent binary histories (figure [TT1 and overall merger rates (figure [T"3t. 
the resultant velocity distributions are very similar. This is a result 
of the ~ t" 1 merger time distribution at late times discussed previ- 
ously. 

The results in figure [14] naturally lead the question of what 
fraction of mergers is consistent with the data to be addressed. 
We wish to assess the fraction of high-mass galactic white dwarfs 
formed by binary mergers (BSE) which we parameterise by 8. This 
results in a fraction 1 — from single star evolution (SSE). To assess 
a value of 9 for a given SSE and BSE Monte Carlo realisation we 
first calculate the galactic formation rate of high-mass WDs from 
SSE and BSE in this realisation, which we denote Tsse and Fbse, 
respectively. Then, for both PG and SDSS we make a copies of 
the BSE objects simulated as observed, and (3 copies of objects 
simulated as observed from SSE. Assuming that equal numbers of 
objects were simulated in both the BSE and SSE realisations, then 
the two simulated samples combined have a galactic BSE fraction 
of 

a _ /3F B se 

PTbse + arssE 

To test whether the data is consistent with t his realisatio n, we 
use the two sample Anderson-Darling statistic dPettitj[T"976h . The 
Anderson-Darling test considers the difference between the sam- 
ples across the entire distribution, and so is more statistically 
powerful that the more commonly used Kolmogorov-Smirnov test 
which depends only on the extremum. The number of simulated 
white dwarfs is always much larger, by at least a factor of ten, than 
the number observed. 

The results for one particular choice of SSE and BSE model 
are shown in figure [T3(aT| In figure [L5"(b)| we show the combined 
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is the case, since when the SDSS objects are plotted in the U — 
V plane (assuming zero radial velocity) 7 of the 9 objects lie in 
the negative U, negative V quadrant. Depending on the unobserved 
radial velocity, many of these could have kinematics consistent with 
the Pleiades and Hyades moving groups. Indeed it has been shown 
that the the white dwarf GD 50 , has a velocity an d cooling age 
consistent with a Pleiades origin l lDobbie et alj|2006h . 

That the data rules out a white dwarf merger origin for the 
majority of high-mass white dwarfs appears more secure, despite 
the apparent consistency of the SDSS sample with the BSE simu- 
lations: The PG sample is entirely consistent with SSE, and nei- 
ther sample contains a high-mass white dwarf travelling at > 
50 kms -1 which would be convincing evidence of a BSE ori- 
gin for some high mass white dwarfs. This is not surprising, since 
the expected number of merger products observed in PG and SDSS 
(^bse in tabled is significantly smaller than the observed number 
of objects. 

We note that a simpler empirical test for the origin of the high- 
mass white dwarfs is suggested by figure [4] The distribution of 
high-mass white dwarfs is consistent with the velocity distribution 
of the intermediate group that displays the kinematics of young ob- 
jects at the 13 per cent level by the Anderson-Darling test. This 
ignores the selection effects which the Monte Carlo simulation ad- 
dresses, but does suggest that the entire combined group of high 
mass white dwarfs is broadly consistent with SSE. 



Figure 14. Predicted distribution of transverse velocities observed in SDSS 
and PG resulting from the merger of CO+CO WDs with 0.95 < Mi + 
M2/M0I.4. Lines are the BSE code with acE^ = 2 (solid blue line, 
model i), the BSE code with oceA = 1 (dotted blue line, model ii) the 
SeBa code using the ■ya common envelope prescription (dashed blue line, 
model iii) and the SeBa code using the aa prescription (dash-dot blue line, 
model iv). Both SeBa models use «ce-^ = 2 and 7 = 1.5. The red line is 
the predicted distribution of transverse velocities resulting from single star 
evolution to a 0.95 < M/Mq < 1.4 white dwarf according to model A in 
table [5] and the black line are the observed distributions. All BSE models 
use a constant SFR and the disk heating relation of model A in tablef5] 



probability that the PG and SDSS samples are consistent with each 
value of 9. In table[5]we summarise the results of this procedure for 
the range of the BSE and SSE models described in tables |3]and|4] 

The results in table [5] show that for the majority of models 
the fiducial numbers of white dwarfs formed via SSE and BSE are 
consistent with the data. The results taken at face value would also 
appear to show that, for most models, at the 1 per cent probability 
level, high-mass white dwarfs must come from a combination of 
single star evolution and mergers of high-mass white dwarfs. This 
appears artificial however: From the right column of figure |15(a)| 
the PG sample is consistent with all SSE, while the SDSS sample 
that has a low probability of arising purely from SSE. 

This is a result of the lack of low velocity (< 14kms -1 ) 
white dwarfs in the SDSS sample. It may be that the lack of low 
velocity white dwarfs in SDSS is a statistical anomaly, since the 
number of objects is small. In theory this would be taken account 
of in the analysis described above, however young stellar objects 
can display prominent substructure in their kinematics as a result 
of moving groups (e.g. lDehnenlll998ri . This would have the result 
of both reducing the effective sample size, and producing a very 
different velocity distribution than the Schwarzschild distribution 
assumed in the SSE Monte Carlo. There are indications that this 



6 SCALE HEIGHTS 

One of the key results of this study is that hot white dwarfs 
of mass > 0.75 Mq had much shorter main sequence lifetimes 
than their lower mass counterparts, and hence their kinematics 
are characteristic of young stars. A direct result of this is that 
these higher mass white dwarfs will have reduced scale height. 
This is vitally important to consider when calculating the forma- 
tion rate as a function of mass us in g local samples suc h as in 
iLiebert. Bergeron & Holberd d2005h or lKepler et"atl d2007l) or pro- 
ducing galactic white dwarf simulations such as Nelemans et alj 

feooih . 

Unfortunately, neither the SDSS or PG sample allow accurate 
direct determination of the scale height of each white dwarf popu- 
lation, particularly the rare and less luminous high-mass groups. In- 
stead, here we list the expected scale height by comparison with the 
SSE models that appear to accurately describe the kinematics. We 
do this to allow simple initial corrections without resorting to the 
simulations of the type performed in this work. The scale height, h, 
was defined through 



v{z) = vq sech 



2 (—) 
\2h) 



(12) 



where u(z) is the stellar number density in terms of the height 
above the plane of the galactic disk, z. The scale height, h, was es- 
timated by constraining equation[T2]to give both the correct overall 
number and central WD density, Vo- We choose this method since 
the most common usage of the scale height is to calculate galactic 
birthrates from local densities. The results are give in table[6] Note 
that the higher mass groups smaller scale height results in a local 
density enhanced by more than a factor of two over the more com- 
mon low-mass group. In particular, the apparent excess of high- 
m ass white dwarfs found in the PG su rvey (discussed in section 6 
of iLiebertB ergeron & Holberg 2005) can be naturally explained 
by their lower scale height, which causes a high abundance in this 
relatively local survey. That the number of high-mass white dwarfs 
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(a) The left column shows the cumulative distribution of transverse veloc- (b) The combined probability that each value 8 is consistant with both the 
ities of high-mass (M > O.95M0) white dwarfs in the SDSS and PG PG and SDSS samples. Calculated by the product of the probabilities in the 
survey. The dashed-dot lines are the predictions of SSE model C and the right column of figure [T"5(a)1 
dashed lines are the predictions of BSE model iv. The right column shows, 
for each fractional galactic formation fraction from BSE, 8, the probability 
that the velocity distribution is consistent with the data using the Anderson- 
Darling statistic for the PG sample, PpQ, and the SDSS sample, Psdss- 
The fiducial 9 is the fiducial predicted galactic fraction from BSE model iii 
compared to SSE model C with 50 per cent binary fraction. 

Figure 15. Plots showing the calculation of the galactic formation fraction of high-mass white dwarfs formed in mergers during binary star evolution in model 
C compared to single star evolution model iii. 
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Table 5. Summary of the results of calculation of the fraction of high-mass white dwarfs formed in mergers compares to single star evolution. The SSE models 
are described in tableland the BSE models are described in table|4] Tbse i s the galactic formation rate (in yr - 1 ) from binary star evolution assuming that the 
merger of two CO white dwarfs with combined mass between 0.95 Mq and 1.4 Mg results in a high-mass white dwarf. Tsse is the galactic formation rate 
from single star evolution. 8 is the galactic fraction of high-mass white dwarfs formed from BSE so that the fiducial value is given by 8^ = -p — rs pp . 

^ BSET-i SSE 

The numbers ATgSE and TVbse are the predicted observed numbers from SSE and BSE evolution respectively in the PG and SDSS samples. P(8fid) is the 
probability that both the PG and SDSS velocity distributions are consistent with 8^ using the Anderson-Darling statistic. 8(P > 0.05) is the range of 8 
values which have a probability of being consistent with the data greater than 1 per cent. The fiducial value of 8 is calculated assuming a 50 per cent binary 
fraction (i.e., two-thirds of all stars formed in binaries). Both SSE and BSE models use the same disk heating model and star formation history: model C of 
table[3]for the constant SFR models, and model D for the exponential. 



is consistent with single star expectations in PG is confirmed by the 
number of expected white dwarfs from single star evolution in table 
® 



7 SUMMARY 

We have analysed the kinematics of young (< 3 x 10 8 years) 
DA white dwarfs from both the PG and SDSS surveys and find 
a strong connection between their mass and kinematics: low-mass 



Table 6. Scale heights, h, defined through equation \l2\ for three different 
mass groups, h is calculated by matching the central density and overall 
number to the simulations described in section l4~2l 



Mow / Mq 


M high / M 


h 1 pc 


0.45 


0.75 


120 


0.75 


0.95 


58 


0.95 


1.40 


54 
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white dwarfs (0.45 M Q < Mi +M 2 < 0.75 M Q ) display the kine- 
matics of old stars, with higher velocity dispersion (~ 46 km s -1 ) 
and asymmetric drift, while higher mass white dwarfs (0.75 Mq ^ 
Mi + Mi < 0.95 Mq) display the kinematics of young stars with 
a velocity dispersion of only ~ 19 km s~ 1 . We have shown in sec- 
tion|4]that this is expected due to the shorter precursor lifetime of 
the more massive progenitors, and that there is agreement both on 
simple analytic grounds (section |4~"0 and more quantitive Monte 
Carlo simulations of the PG and SDSS samples (section |4~2l l. 

A further key conclusion is that the white dwarf scale height 
and its variation with age and mass is vitally important to consider 
when calculating birth rates based on local samples (section|6]l. 

In addition, we have separately analysed the highest mass 
white dwarfs (M > O.95M0, section [5), since it has been sug- 
gested that many of these formed as a result of the merger of two 
lower mass CO white dwarfs. We find at present a discrepancy in 
the SDSS velocity distribution where no high-mass white dwarfs 
with transverse velocity less than 14 kms -1 is detected. This re- 
sults in a velocity distribution that within our statistical framework 
is inconsistent with purely single star evolution. We argue this is 
likely to an anomaly, either be a statistical, or a result of a number 
of these white dwarfs being members of moving groups. We find 
that, even under the most optimistic binary evolution models, we 
would only expect to find 3 white dwarfs formed via white dwarf 
binary mergers and that the apparent excess of high mass white 
dwarfs found in PG is caused by their reduced scale height. In ad- 
dition, we note the kinematic 'smoking gun' of some fraction of 
high-mass white dwarfs coming from binary evolution would be 
high-mass white dwarfs traveling at > 50 kms -1 , of which none 
are found in PG or SDSS. 
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a nuisance parameter, is an ellipsoid in the sky plane 

p(v h v b ) = C" exp [ - a(vt - v t ) 2 - /3(v b - v b ) 2 

- 2j(v e - v e )(v b - Vb)] , (A2) 

where v~i and v~i are the components of Vb in the directions of I and 
b (which can be obtained via (v~i, v~b, vj) = R ■ Vb) and a, /3, 7, 
and C' are given by 



a = A22 — A12/A11 , 
P = A 33 - A? 3 /An , 
7 = A23 - A12A13/A11 , 



C' = \f a/3 — 7 2 /7r . 



(A3) 
(A4) 
(A5) 
(A6) 



For each object we have measurements of vi and v b , together 
with an associated velocity error a. Integrating over the 'true' vi 
and Vb gives the log likelihood used in equation[9]as 

log Ci{v?",v?")= log J dVf(V)P(V\v? s ,vZ hB ,a) 
= logC" 



(a + S)(P + 6) -7 2 
[(Av 2 b + Av!)(a/3 - 7 2 )+ 
S(pAv 2 b + aAvj + 2-yAvtAvb)] , (A7) 



where 



S = l/2cr 2 , 

Aobs — 
Vi = v e - Vi , 

Aobs — 
Vb = V b — Vb , 

C" = C 



= 5 



a/3 — 7 2 



(A8) 
(A9) 
(A10) 

(All) 
(A 12) 



Note that for small error, S — > 00, and equation I A7l reduces to the 
log of equation lA2l as expected. 



APPENDIX A: LIKELIHOODS 

Here we give our expressions for the proper motion like- 
lihoods of an individual object. These largely follow 
iRatnatunga. Bahcall & Casertanol ( I1989h . modified to include 
errors in proper motion. We ignore errors in sky position (£, b), 
which are small. 

Assuming a Schwarzschild distribution function, then, in co- 
ordinates aligned with the principle axes of the velocity ellipsoid, 

f(V) = exp (-(V - v ) T • r • (V - Vo)) , 

(Al) 

where T = diag(l/2<ri , 1 /2<T2, 1/20-3) and Vb is the mean veloc- 
ity. Ignoring errors in distance, we then rotate to axes aligned with 
the sky plane, and integrate over the unobserved radial velocity, 
which, in this case, is a nuisance parameter. 

We define, A, to be the dispersion tensor rotated into the co- 
ordinate system, (I, b, d), aligned with the sky plane. This will be 
given by A = R ■ T , where R is a rotation matrix (given ex plic- 
itly as equation A4 in lRatnatunga. Bahcall & Casertandl989h . The 
probability distribution, after integrating over the radial velocity as 
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